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ABSTRACT 

Several recent observations suggest that gas-poor (dissipationless) mergers of elliptical galaxies 
contribute significantly to the build-up of the massive end of the red sequence at z < 1. We perform a 
series of major merger simulations to investigate the spatial and velocity structure of the remnants of 
such mergers. Regardless of orbital energy or angular momentum, we find that the stellar remnants lie 
on the fundamental plane defined by their progenitors, a result of virial equilibrium with a small tilt 
due to an increasing central dark matter fraction. However, the locations of merger remnants in the 
projections of the fundamental plane - the Faber-Jackson and R e — M» relations - depend strongly 
on the merger orbit, and the relations steepen significantly from the canonical scalings (L oc o\ and 
R e oc M®- 6 ) for mergers on radial orbits. This steepening arises because stellar bulges on orbits with 
lower angular momentum lose less energy via dynamical friction on the dark matter halos than do 
bulges on orbits with substantial angular momentum. This results in a less tightly bound remnant 
bulge with a smaller velocity dispersion and a larger effective radius. Our results imply that the 
projections of the fundamental plane - but not necessarily the plane itself - provide a powerful way 
of investigating the assembly history of massive elliptical galaxies, including the brightest cluster 
galaxies at or near the centers of galaxy clusters. We argue that most massive ellipticals are formed 
by anisotropic merging and that their fundamental plane projections should thus differ noticeably from 
those of lower mass ellipticals even though they should lie on the same fundamental plane. Current 
observations are consistent with this conclusion. The steepening in the L — a e relation for luminous 
ellipticals may also be reflected in a corresponding steepening in the Mbh — c e relation for massive 
black holes. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation 
— galaxies: fundamental parameters 



1. INTRODUCTION 

Several independent lines of evidence point to the 
importance of gas-poor (dissipationless) mergers in 
the assembly of massive elliptical galaxies. For ex- 
ample, the decreasing rotational support and transi- 
tion from disky to boxy isophotal shapes with in- 
creasing stellar mass in elliptical galaxies suggests 
an increasing fraction of dissipationless mergers in 
the g rowth of the most massive elliptical galaxies 
/e.g.. IBenHer et al.l iKormendv fc Bended EH 

TB1997H N 



l)van Dokkum et all 19991: iTran et al.ll2005t) as well as i n 



iFaber et alJll997tlNaab et alJl2006T). The lack of correla- 
tion between me tallicitv llGallazzi et al"l l2005') or metal- 
licity gradients ijCarollo et al.1 119931) and galaxy mass 
for massive ellipticals, while both are positively corre- 
lated with galaxy mass for lower-mass ellipticals, pro- 
vides additional support for this picture. There is 
also indirect eviden ce from the evolution of galaxy lu- 
minosity func tions (|Bell et all 120041: IFaber et all 120051 
lBlantonl2005l) and semi-analytic models of galaxy forma- 
tion lIKauffmann fc Haehneltil200Ct IKhochfar fc Burkertl 
2003t|De Lucia et aLll2006(l that a significant fraction of 
low-redshift ellipticals were assembled via dissipationless 
mergers. Moreover, there is direct observational evidence 
for a substantial number of "red" mergers between early- 
type galaxies in an intermediate redshift galaxy cluster 
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more local galaxies (B ell et al .1 120051 Ivan Dokkumll27) 05) . 

In spite of the importance of gas-poor mergers, however, 
the properties of these mergers and their remnants are 
not fully understood. 

If elliptical galaxies gain substantial mass via dissi- 
pationless merging at z < 1, this process could have 
a significant impact on their properties. The purpose 
of this paper is to study how such mergers affect the 
global dynamical and kinematic structure of the rem- 
nants, as ellipticals show a very tight correlation in the 
space spanned by their half-light radii (i? e ), central ve- 
locity dispersions (<r e ), and central surfa ce brightnesses 
(I F ) k nown as the fundamental plane ({Dressier et alJ 
119871 iDiorgovski fc Davisl 119871) . Relationships be- 
tween two of the fundame ntal plane variables, su ch as 
the Faber-Jackson relation (|Faber fc Jackson|ll976|) and 
the stellar mass-size relation, may also serve as im- 
portant constraints for galaxy formation models (e.g., 
IMcIntosh et aDl200l . 

This is particularly true for the most luminous ellipti- 
cal galaxies, including brightest cluster galaxies (BCGs) 
that populate the exponential tails of the galaxy lumi- 
nosity, mass, and velocity dispersion functions and con- 
tain the most massive black holes in the universe. These 
galaxies are likely to have undergone multiple generations 
of gas-poor mergers, as repeated merging during cluster 
formatio n is a well-motivated mechan ism for BCG for- 
mation l|Merritrill985l lDubinskl H9981. In this model, 
anisotropic infall along cosmological filaments introduces 
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a pre ferential direction in the merging process ijWesti 
1994), in contrast to mergers of galaxies in less ovcrdensc 
environments. As discussed below, we suspect that this 
special assembly history has implications for the scaling 
relations of massive ellipticals: the preferentially radial 
orbits of mergers that form the most luminous ellipticals 
may lead to different Faber- Jackson and stellar mass-size 
relations for these galaxies than for normal ellipticals. 

There are, in fact, tantalizing observational hints of 
deviations from the canonical elliptical galaxy scaling re- 
lations for BCGs and othe r lumi nous ellipticals. For ex- 
ample, lOegerle fc Hoessel 1)1991)) find that while BCGs 
lie on the same fundamental plane as lower mass ellipti- 
cals, their central velocity dispersions change very little 
with increasing luminosity and their effective radii in- 
crease steeply with luminosity, significantly more so than 
for normal ellipticals. A similar effect is present in the 
ea rly-type galaxy sa mple of the Nuker team (see fig. 4 

of iFaber et alJ 119971 and Lauer et al. 2006, in prep.) 

as well as in th e iPostman fc Lau erl l)1995l) BCG sample 
l)Grahamlll996j) . while iWvithd 1)2006) has claimed to de- 
tect an analogous change in the Mbh — c e relation for 
the most massive black holes. 

In this paper, we present the results of a large num- 
ber of simulations of major mergers of elliptical galaxies. 
Our galaxy models contain dark matter haloes and stel- 
lar bulges, both simulated with sufficient force and mass 
resolution to reliably investigate the remnants down to 
scales of ~ 0.1i? e (over 10 6 total particles per simula- 
tion). We consider a wide variety of energies and angu- 
lar momenta for the merger orbits and quantify the ef- 
fects of orbital parameters on the merger remnants. This 
work is a significant extensio n of our earlier related study 
l)Bovlan-Kolchin et al.ll2005l hereafter BMQ), where we 
considered only a few possible orbits for the merging 
galaxies. The simulations presented in this paper allow 
us to study in detail the effects of dissipationless merg- 
ing and merger orbits on the fundamental plane and its 
projections. Section |2 describes the initial galaxy mod- 
els and simulation parameters. Section|21contains the re- 
sults of our simulations, including the locations of merger 
remnants in the fundamental plane fSection l3.1)l and its 
projections (Section EOfl . In Section^] we discuss our re- 
sults and their implications for the assembly of massive 
elliptical galaxies and their central black holes. 

2. METHODS 

The galaxy models used in our simulation consist 
of dark matter haloes and stellar spheroids and are 
fully described in BMQ; we review the essential fea- 
tures here. The in itial dark matter haloes have the 
iNavarro et alJ l)1997)) density profiles with a virial mass 
M200 = 10 12 Mq and a concentration c = 10, resulting in 
a virial radius of 162.6 kpc and sc ale radius of 16.26 kpc. 
The initial stellar bulges have the lHernauisti 1)1990)) den- 
sity profiles with a total stellar mass M* = 5x 10 10 Mq, 
resulting in a stellar baryon fraction of /{, = 0.05. The 
additional free parameter in the stellar profile, the scale 
radius a, is set by the R P — relation m easured for 
SDSS elliptical galaxies bv lShen et all 1)200.11) : 

/ ji ,r \ 0.56 

^ = 4 - 16 (lo^J ^ « 



where R e — 1.8153 a for the Hernquist profile; this gives 
R e = 2.82 kpc for our initial models. This relation is 
appropriate for low-redshift elliptical galaxies. Ideally, 
we would use the R e — M* relation for z w 1 rather 
than z » 0. The data that are available at higher red- 
shift indicate there has not been strong e volution in the 
R e — relation from z = 1 to z = (Tru iillo et alJ 
I2004U iMcIntosh et alJl2005|) . We include the effects of 
b aryonic dissipation via the adiabatic contraction model 
of iBlumenthal et alJ 1^)86). Although the applicability 
of this model to elliptical galaxies is not fully obvious, 
recent hydrodynam ic simulations sugge st that it is sur- 
prisingly accurate l)Gnedin et al.ll2004)) . The resulting 
(circular) aperture velocity dispersion within R e for the 
stellar component in our models is a e — 151 km/s. We 
also perform a number of unequal mass mergers, with 
mass ratio 0.33:1, where the less massive galaxy also has 
fb = 0.05 and a scale radius set by equation QJ. Al- 
though we present the parameters of our initial condi- 
tions in physical units, the simulations reported here can 
be readily scaled to other stellar and dark matter masses 
(see associated disc ussion in Sect ion 13.211 . 

We use GADGET- 1 l|Springel et"ail20 01) to perform the 
merger simulations. All runs use a force softening of 
0.3 kpc (= 0.106i? e ) and equal-mass particles for the 
stellar and dark matter components: N* = 2.5 x 10 4 , 
Nbm — 5 x 10 5 for each galaxy model. The initial par- 
ticle positions of each component are sampled from the 
component's spherically symmetric density profile. To 
ensure that the models are in equilibrium, we compute 
separate distribution functions for the halo and bulge and 
then use these distribution functions to initialize the par- 
ticle velocities. Our initial ga laxy models do not include 
any velocity anisotropy. As INipoti et alJ l)2002f) show, 
however, the inclusion of a modest amount of velocity 
anisotropy in models of merging galaxies does not sub- 
stantially change the properties of the remnants, while 
adding significant velocity anis otropy mak e s the initial 
model unstable. Additionally, iCiotti et alJ l)1996)) show 
that velocity anisotropy cannot be the primary origin for 
the tilt in the fundamental plane. 

iV-body models set up with this procedure are in 
gener al extremely stable ov er many dynamical times 
Ce.g.. IKazantzidis et alJl2004j) . which we have explicitly 
checked for our models. An additional source of spurious 
evolution inherent to ./V-body simulations is two-body 
relaxation, a consequence of using simulation particles 
that are many orders of magnitude more massive than 
the stars or dark matter particles they represent. We 
have performed test runs to ensure that two-body relax- 
ation docs not affect our results for the duration of the 
simulations (~ 6 Gyr) for our choices of particle number 
and force resolution (see BMQ for detailed resolution and 
stability studies). 

The centers of the galaxies in the merger simula- 
tions are initially separated by the combined virial 
radii of the two haloes and the orbital parameters 
are defined in the standard two-body point mass ap- 
proximatio n (Bin nev"fc Tremaind 119871 Appendix ID; 
iKhochfar fc Burkertll2006)) . In this limit, orbits can be 
characterized by two quantities: energy E and angular 
momentum L, or equivalently, pericentric distance r p and 
eccentricity e (or circularity e = vl — e 2 for runs with 
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TABLE 1 

Summary of Merger Simulations 



Run a 


Tl IT h, 

Mass b 


Orbit c 


~\T d 
» orb 




f 


PI 


1*1 


l arabolic 




n 
u 




P2 


1*1 


Parabolic 


9^n 

ZJU 


2 29 




P3 


1:1 


Parabolic 


250 


12.5 


- 


Pi 


1:1 


Parabolic 


250 


37.8 




P5 


0.33:1 


Parabolic 


222 







P6 


0.33:1 


Parabolic 


222 


2.31 




P7 


0.33:1 


Parabolic 


222 


12.5 




P8 


0.33:1 


Parabolic 


222 


20.5 




Bl 


1:1 


Bound 


200 





0.00 


B2 


1:1 


Bound 


200 


2.29 


0.10 


B3 


1:1 


Bound 


200 


20.8 


0.30 


B4 


1:1 


Bound 


200 


37.8 


0.40 


Bo 


1:1 


Bound 


200 


57.1 


0.48 


B6 


1:1 


Bound 


200 


78.4 


0.56 



a Name of run 
b Mass ratio of progenitors 
c Energy of run (either parabolic or bound) 
d Initial orbital velocity (km s~ 1 ) for a galaxy with 
M dm = 10 12 M© and M,=5x 1O 1O M 
c Pericentric distance of orbit (kpc) 
'Circularity of orbit (undefined for parabolic orbits) 



E < 0) . The orbital parameters for all of our production 
runs are listed in Tabled Since our goal is to understand 
the effects of different orbits on the properties of the 
merger remnants, we consider both parabolic encounters 
and bound orbits, and for each energy we perform sim- 
ulations over a range of orbital angular momenta. The 
energy of the bound orbit s is taken to co rrespond to the 
most probable energy from lBensonl 1)20051) , which is based 
on a study of dark matter halo mergers averaged over all 
cosmological environments. We use this value of E only 
as a rough guide and note that existing statistical stud- 
ies of dark matter halo mergers are dominated by low 
mass objects and therefore do not necessarily represent 
the orbits responsible for the assembly of massive ellip- 
tical galaxies in cluster environments. Larger samples of 
haloes should soon provide robust statistical information 
about the environmental dependence of halo mergers. In 
particular, as we argue later in this paper, major merg- 
ers leading to the formation of central massive galaxies 
in clusters may follow preferentially radial orbits. 

The non-spherical nature of the luminous portion of 
merger remnants means that quantities such as R e and 
a e depend on the orientation at which the objects are 
viewed. An immediate consequence is that an intrinsic 
scatter in the fundamental plane projections will exist 
solely due to a variation in viewing angle. In order to 
compute R e and a e that can be compared with obser- 
vations, and to quantify the effects of viewing angle, we 
observe the remnant in each simulation from I0 4 different 
lines of sight chosen randomly from points on a sphere. 
The effective radius is determined by the radius of the 
circle containing half of the projected stellar mass and 
the mass-weighted aperture velocity dispersion is com- 
puted within one R e . 

3. RESULTS 
3.1. The Fundamental Plane 



The fundamental plane is typically expressed in the 
coordinate system defined by the observables (x, y, z) = 
(logcr e ,log/e,logi? e ) as 



R e cx a a A~ 



z — ax — by + zq 



(2) 



where a represents the angle of the plane from the ver- 
tical axis z — logi? e , and b/a represents the amount of 
rotation about the z-axis. If we write 



M*(< R e ) = ci I e R 2 e and M dyn (< R e ) = c 2 



G 



( 3 ) 

where c± and c 2 are structure coefficients, then the fun- 
damental plane observables are related by 



R e oc — 

Cl 



C2 /M, 



i dyn 

A/, 



(4) 



ijBender et al.ll 992). Therefore, if elliptical galaxies form 
a homologous sequence (if c\ and c 2 are constant) with 
light tracing dynamical mass and a constant stellar mass- 
to-light ratio, they should lie on the virial plane given by 
a = 2, b = 1. However, the observed fundamental plane 
differs from the virial plane. For example , a = 1.24±0.07. 
b = 0.82 ±0.02 in the Gunn r-band Iporgensen et alJ 
Il99l and a = 1.53 ±0.08, b = 0.79 ±0.03 in the K-band 
ijPahre. Diorgovski. &: de Carvalh"olll998b|) . Results for 
all SDSS wavebands are similar t o the K-band fundamen- 
tal plane ijBernardi et al.l|20 03a'). This difference repre- 
sents a "tilt" in the observed fundamental plane from 
the virial plane: either c^/ci, -Md yn /Af*, or M*/L varies 
systematically with galaxy mass. 

We are able to investigate the first two possibil- 
ities with our simulations, but we cannot directly 
test the third. In what follows, we assume a con- 
stant stellar mass-to-light ratio for the galaxies in our 
simulations. We can therefore interpret our results 
as testing the amount of the fundamental plane tilt 
that is due to either non-homology or central dark 
matter fraction. A number of authors have argued 
that stellar population variations with galaxy mass 
are responsible for the fundamental plane tilt (e.g. 



Gerh ard et a 
20061) . 



As 



1120011: FTreu et all 120051: iCappellari et al 
Pahre. de Carvalho. fc Dioreovskil l)1998 



point out, however, the existence of a significant tilt in 
the near-infrared, where stellar population effects should 
be less pronounced than in the optical, indicates that 
stellar populations alone are not the cause of the fun- 
damental plane tilt. Furthermore, K-band M^/L ratios 
do not vary strongly with the age of the stellar popula- 
tion for relatively old stellar populations that are ap- 
propriate for the galaxies we are simulating. Analy- 
ses of SDSS early-type gal axies l|Kauffmann et alJl2003t 
iPadmanabhan et alJ 12004^) also find that stellar mass- 
to-light ratios are relatively constant with mass. We 
therefore believe that assuming a constant stellar mass- 
to-light ratio is a reasonable approximation. 

Fig. ^ shows the locations of all of our equal-mass 
merger remnants in the edge-on projection of each of 
the three planes discussed above - virial, r-band, and 
K-band. We find that the 0.33:1 merger remnants show 
all of the same trends as the 1:1 remnants and thus the 
former are not plotted here. The 1 — a errors of the fun- 
damental plane in the direction of R e are indicated by 
the dotted lines in Fig. ^ Each remnant is marked with 
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5.0 5.2 5.4 3.1 3.2 3.3 3.4 3.5 3.6 3.8 3.9 4.0 4.1 4.2 4.3 

2-0 log l0 (O-1.0 log 10 (I.) 1.24 log 10 (ff.)-0.82 lo gl0 (I e ) 1.53 log 10 (a e )-0.79 log 10 (I e ) 

Fig. 1. — Locations of our equal-mass merger remnants in the plane in (log a e , log I e , log R e )-space that corresponds to the edge-on 
projection of the v irial plane (left), r-band fundamental plane (middle; IJorgensen et all 119961) . and K-band fundamental plane (right; 
IPahre ct al. 1998tj). The solid line in each panel indicates a slope of 1 for the edge-on view; the dotted lines in the middle and right panels 
indicate the observed l — cr scatter in the direction of R e . For each panel, each simulation (corresponding to a particular energy and angular 
momentum) is plotted with a different colour. Each point of a given colour corresponds to one of 10 4 viewing angles, while the square 
marks the "most probable" remnant. A single re-merger simulation, using the remnants of runs P3 and B6 as initial galaxies and with a 
parabolic orbit with r p = 21 kpc, is plotted in black in the upper right corner as well. The initial galaxy model is a single point (large 
triangle). The effective radii are measured in kpc, the velocity dispersions in km s _1 , and surface brightnesses in 10 10 Mq kpc - 2 for all 
panels. 
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Fig. 2. — Locations of the merger remnants from Fig. Ill in (log M^ yn , log M», log i?, e )-space: M^ yn oc M^^R". Both M<j yn and M» are 
computed within a sphere of radius R e . Left: The v = plane, which is a good approximation to the r-band plane (fig. lb). The dotted 
lines correspond to M^ yn oc M t +fl with fi = 0.12 (lower) and fi = 0.25 (upper); these bracket the results of our simulations. Right: The 
fi = plane, which fits the observed K-band fundamental plane well. We plot the parameterization of the tilt in the fundamental plane 
that best fits our simulations: M^ yn /M f <x R®' 21 (dotted line). The effective radii are measured in kpc, while the stellar and dynamical 
masses are measured in Mq. 



a different colour, and each point represents one random 
viewing angle of the remnant. The initial galaxy is spher- 
ical by construction, so it is located at a single point in 
the fundamental plane (large triangle; lower- left corner). 
The stellar mass-to-light ratio has been chosen so that 
the initial model lies on the observed fundamental plane 
(this corresponds to K-band and r-band mass-to-light ra- 
tios of 1.0 and 3.6, respectively); the same stellar mass- 



to-light ratios are used to place the merger remnants on 
the fundamental plane. 

The left panel of Fig. ^ shows that all of the remnants 
lie significantly off the simple virial theorem expectation. 
Strikingly, however, all of the simulated observations of 
all of the merger remnants lie well within the l — cr errors 
of the fundamental plane for both the optical and near- 
infrared planes. These results establish that the rem- 
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nants of major mergers of spherical elliptical galaxies lie 
on the fundamental plane defined by their progenitors for 
the wide range of orbital energies and angular momenta 
covered by our simulations. We note, however, that our 
dynamic range in R e and M* is rather modest, and all 
of our simulations begin at the same point in the funda- 
mental plane. Thus our simulations cover only a limited 
part of the observed plane. 

Since massive elliptical galaxies may have undergone 
multiple dissipationless mergers during their assembly 
history, we also performed a simulation in which the two 
progenitors were randomly chosen to be the remnants 
of runs P3 and B6 and the orbit was parabolic with 
r p = 21 kpc. We find this re-merger remnant lies on the 
same fundamental plane as the initial galaxy model and 
first-generation merger remnants (Fig. ^ black region in 
upper right corner). This single re-merger simulation is 
certainly not definitive, but it does support the hypoth- 
esis that the presence of elliptical galaxies on the funda- 
mental plane is consistent with multiple generations of 
dissipationless merging. 

While the observed fundamental plane is typically 
plotted in the (x,y,z) = (logfj e ,log/ e ,logi? e ) space 
as shown in Fig. 1, we find it useful to interpret the 
plane theoretically in the transformed coordinate system 
(x',y',z') — (logMdy n) log -M*, log R e ), in which a plane 
can be written as 

M dyn cx M^R" , or z' = ± [x' - (1 + + z' , 

(5) 

where both Md yn and M* refer to the masses enclosed 
within R e . The parameters [i and v quantify how the 
ratio of dark matter mass to stellar mass within R e varies 
with increasing galaxy mass and size. 

Comparing equations @ and (JSJ), and using the defi- 
nitions in equation gives 



The virial plane of a = 2, b = 1 corresponds to fj, = 
0, v = (i.e. Afdyn oc M*); it is therefore conveniently 
parallel to axis z' — \ogR e in the Afdyn — M* — Re 
space. The observed r-band and K-band fundamen- 
tal planes have (/i, v) = (0.323, —0.0323) and (/i, v) — 
(0.0327, 0.242), res pectively. It is interesti ng to note that 
the r-band plane of Jor gensen et alJ (1996) has v sa 0, i.e. 
Mdyn/M* oc M*, while the K-band fundamental plane 
has n w 0, i.e. M dy „/M* oc R v e . 

Fig. 2 shows the locations of our merger remnants in 
the edge-on view of the v = (left panel; appropriate for 
r-band) and [i = (right panel; appropriate for K-band) 
planes. The remnants appear more consistent with the 
K-band plane than the r-band plane. Although Fig. 2 
and Fig. lbc display the same information in different 
coordinate systems, Fig. 2 illustrates clearly how the dy- 
namical mass within R e increases more than the stellar 
mass within the same radius, meaning that the inner 
regions of our merger remnants are more dark matter- 
dominated than their progenitors. This increasing dark 
matter fraction with increasing galaxy mass is the pri- 
mary cause for the tilt of the observed fundamental plane 
from the virial plane in our gas-free mergers. 

Since systematic non-homology and varying cen- 
tral dark matter fraction in elliptical galaxies are 



often considered as possible ori gins for the tilt 
of the fundame ntal plane (e.g., iFaber et alJ 119871 
I Ciotti et al.lll99ll IGraham fc Collesslll997t IPahre et all 
| 1998al iGerhard et alJ 120011: iPadmanabhan et all 120041 
Trui illo et al.l2004aj) . it is useful to assess the importance 
of non-homology in the context of the two-component 
dark matter and stellar models considered in this paper. 
We find that both c\ and C2 [as defined in equation 
are constant for our merger remnants 3 : c\ = 2.70 with a 
maximum deviation of < 4% and C2 = 2.90 with a maxi- 
mum deviation of < 3%. This implies that the final sys- 
tems (bulge+halo) are in virial equilibrium and are dy- 
namically homologous (because ci is constant) and that 
the luminous portions of our remnants are structurally 
homologous (because c\ is constant). The total system is 
structurally non-homologous, however, due to a change 
in the dark matter halo relative to the bulge. The inner 
part of the halo, which contains the bulge, does not grow 
in radius as much as the bulge does during the course 
of each merger, while the outer part of the halo expands 
significantly. This results in an increased central dark 
matter fraction for the remnants, as shown in Fig. 2. 
The amount of tilt in our simulations is very similar to 
the observed fundamental plane tilt, implying that there 
is no need to invoke stellar mass-to-light variations to 
account for the tilt in our simulations. 

As the bulges themselves are structurally homologous 
(ci is constant and the density profiles of the remnant 
bulges are well-fit by a Hernquist profile), we refer to 
a varying dark matter fraction rather than a structural 
non-homology as th e origin of the tilt. We note that 
IGraham et al.1 l)1996j) find Sersic r 1 /™ models provide bet- 
ter fits to BCG surface brightness profiles than the r 1 / 4 
law and that n correlates positively with mass, indicat- 
ing a real non-homology in the elliptical galaxy popu- 
lation that is not ca ptured by our simulations (but see 
iGonzalez et al.ll2005l for an argument that two separate 
r 1 / 4 profiles fit BCGs significantly better than Sersic pro- 
files). 

A further point of interest is that the scatter due to 
viewing angle variation of the non-spherical remnants in 
our simulations (represented by the spread within each 
colour in Fig. ^) is roughly parallel to the observed fun- 
damental plane. This trend is a consequence of the 
virial equilibrium of individual galaxies: at fixed stel- 
lar mass, (igi? e ~ i?g~ - 33 , so R e oc er e ' 2 ~ 3 * ) . Intuitively, 
a non-spherical remnant maintains a small line-of-sight 
velocity dispersion by extending in the perpendicular 
direction. This shows that while viewing angle varia- 
tion introduces some scatter into the fundamental plane 
l)Gonzalez-Garcia fc van Albadal200"3^ . this scatter tends 
to be along the fundamental plane itself. In addition, 
the scatter due to viewing angle variation does not in- 
crease between the first and second generation mergers, 
further indicating that the fundamental plane scatter is 
relatively insensitive to dissipationless merger history. 

Similar conclusions about mergers and remergers pre- 
serving the fundamental plane hav e been reached in var- 
ious settings by previous authors l|Canelato et al.lli995l: 

3 Note that we compute the surface mass density, while it is the 
surface luminosity density that is measured observationally. The 
two are simply related by the stellar mass-to-light ratio, which we 
assume is constant. 
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Nipoti et alJl2003t iGonzalez-Garria & van Albad3l200a 
Robertson et alJ 120051) . The previous rem erger simula- 
tions m ost relevant to this work ar e those of lNipoti et al.l 
( 2003) and lRobertson et all <|2005j) . They show that both 
remergers of two-component galaxies similar to ours and 
remergers of triaxial remnants of gaseous disk galaxy 
mergers lie on the (K-band) fundamental plane. These 
results suggest that the spherical symmetry typically as- 
sumed in initial galaxy models of merger simulations, 
albeit unrealistic, does not significantly affect the results 
on the preservation of the fundamental plane. Explor- 
ing the evolution of the fundamental plane under dissi- 
pationless merging for general triaxial initial conditions 
would be useful; however, we believe that a brute-forced 
sampling of the large parameter space associated with 
merging triaxial galaxies is highly inefficient. We will 
discuss a more economical procedure in Section 4 that 
uses the results of cosmological simulations to motivate 
realistic merging sequences. 

3.2. Projections of the Fundamental Plane 

The fact that our merger remnants are consistent with 
the fundamental plane does not necessarily imply that 
the properties of the remnant are insensitive to the 
merger orbit, nor does it imply that the remnants nec- 
essarily follow the observed Faber- Jackson or R e — A/* 
relations: requiring a galaxy of a given mass to lie on the 
mean of both fundamental plane projections picks out a 
single point in R e — a e ~I e space, whereas the fundamen- 
tal plane at fixed mass forms a line in the same space (an 
R e — a e correlation). 

In order to quantify the effects of merger orbits on the 
projections of the fundamental plane, we compute the 
slopes in the R e oc M" and M* cx erf relations for each 
of our 1:1 merger simulations. Fig. [3] shows the depen- 
dence of the slopes a (left) and (3 (right) as a function of 
the pericentric distance of the orbit 4 for both the bound 
orbits (diamonds) and parabolic orbits (squares). The 
slopes have a strong trend with pericentric distance: in- 
creasing r p corresponds to decreasing a and decreasing 
0. These trends were already noted in BMQ based on 
a few orbits, but the comprehensive set of simulations 
with different L and E presented in this paper provide a 
much more quantitative understanding of the impact of 
the merger orbit on the fundamental plane projections. 
The slopes for the re-merger remnant, marked with an 
X in Fig. |3 agree well with those for the first genera- 
tion mergers. Although Fig. shows only the results of 
1:1 merger simulations, we find similar trends in 0.33:1 
mergers. These results suggest that the angular momen- 
tum and energy of the orbit have the biggest impact on 
the fundamental plane projections, with the mass ratio 
of the merging galaxies and the internal structure of the 
initial bulge/halo system being of secondary importance. 

For comparison, the observed fundamental plane 
projections are R e cx £ - 56 -°- 63 anc j L cx c 4-414 
jjPahre et al1ll998bt iBernardi et al]l2003bft . although as 
noted in Section 1, these scalings are probably not appli- 
cable to the most luminous ellipticals. Our high angu- 
lar momentum simulations have scaling relations in good 
agreement with both of the observed projections. The 

4 At fixed energy, varying the angular momentum is equivalent 
to changing the pericentric distance. 



low angular momentum runs, however, deviate signifi- 
cantly from these canonical scalings (similar deviations 
for low angular m omentum r uns were also reporte d in 
iNipoti et alll2003l BMQ, and lRobertson et al.ll20(]l . 

The combination of Figs. Q an d El reveals that while 
all of the remnants lie on the fundamental plane, the 
remnant properties are nonetheless strongly affected by 
merger orbits. Differences in merger orbits can lead to 
significant differences in the projected scaling relations, 
but variations in the R e — M* relation are compensated 
by corresponding variations in the M* — o~ e relation, 
maintaining the fundamental plane because of virial equi- 
librium of the bulge-halo system. 

It is important to point out that although we have 
quoted our results in terms of galaxies with specific 
masses and radii, they can be easily rescaled to other 
masses and radii. Such a rescaling will affect the am- 
plitude of the relations we are studying but not the 
slopes. For example, if we increase all masses by a 
factor of 10 and all radii by a factor of 10 1 / 3 (to pre- 
serve virial scalings), the velocity dispersion will scale as 
^W/R = 10 1 / 3 . This moves our initial conditions off 
of the R e - M* relation, since R e = 2.83 * 10 1 / 3 = 6.1 
kpc rather than the 10.2 kpc given by the Shen et al. 
relation. This difference is comp arable to the ob served 
scatter in the relation, a\ nr w 0.3 l)Shen et all2 003). The 
velocity dispersion will also differ from that predicted by 
the Faber- Jackson relation: er e = 151 * 10 1 / 3 = 325 km 
s _1 rather than 151 * 10 1 / 4 = 269 km s _1 . Thus mergers 
of galaxies at these scales have remnants with identical 
R e — M» and M* — a e slopes as those given above but 
with different amplitudes than the means of the observed 
relations. The slopes of the relations quoted here are 
therefore quite robust. 

Qualitatively, the dramatic steepening of the R e — M* 
and M* — cr e relations with decreasing orbital angular 
momentum in Fig. [3] can be understood by considering 
the role of dynamical friction and energy transfer during 
a merger. At a given orbital energy, stellar bulges with 
substantial angular momentum need to lose more angular 
momentum (and accordingly, more energy) via dynam- 
ical friction on the background dark matter halo than 
do bulges on low angular momentum orbits. This leads 
to a more tightly bound remnant bulge with a larger 
a (smaller (3) and smaller R e (smaller a), as is found 
in the simulations. The same effect also explains why 
orbits with more initial orbital energy (our parabolic or- 
bits) tend to lead to less bound remnants with smaller a 
(larger (3) and larger R e (larger a). 

Quantitatively, we can understand the essence of the 
physics involved in producing the trends seen in Fig.|3|by 
examining the energy equation for the stellar bulges. The 
energy conservation equation for two identical bulges, 
each with initial stellar mass Af*^ and effective radius 
Ri, merging to form a final bulge of mass Af*j and ef- 
fective radius Rf can be written as 



2/i 



Mi 



n- 



M 2 ■ 

*,2 



(7) 



R f Ri 2i?j 

The parameter / depends on the structure of the bulge 
and surrounding dark matter halo: 

Rp. 2 J r W 
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Fig. 3. — Logarithmic slopes of the mass-size relation (R e oc M"; left) and Faber-Jackson (M* oc erf?; right) relations as a function of 
orbital pericentric distance for each of the equal-mass merger simulations listed in Table HI The two orbital energies are shown as diamonds 
(bound) and squares (parabolic), while the re- merger run is marked with an X symbol. At a given T" per ;, the bound orbits produce remnants 
with a slightly lower a and /3 than parabolic orbits. The parabolic head-on orbit has /3 ~ 28, so it is not plotted here. For orbits with 
significant orbital angular momentum (large r per i), our merger simulations reproduce the observed values of a ~ 0.6 and /3 ~ 4 but we 
predict a sharp increase in a and (3 for more radial mergers. The results for the re-merger run agree well with those for the first generation 
mergers, indicating that the spherical and isotropic galaxy models used in the first generation mergers do not strongly bias our results. 



where M^ yn (r) — M*(r) + Md m (r). The parameter 
V = /orb + ft measures the orbital energy of the bulge- 
bulge system when the bulges "touch." r\ can be de- 
composed into two separate parts, with one contribution 
from the initial orbital energy at large radii (/ or b) an d the 
other from the subsequent energy transfer between the 
stellar bulges and dark matter haloes during the merger 
(ft); see BMQ for more details. The values of / and rj 
for our equal-mass merger simulations are listed in Ta- 
ble|21 These results show that, as argued in the previous 
paragraph, at fixed initial orbital energy (fixed forb), in- 
creasing the orbital angular momentum results in a larger 
value of rj, which is due to more energy transfer from the 
bulges to the haloes (a larger positive value of f t ). This 
in turn leads to a more tightly bound remnant bulge 
with a larger a and smaller R e (Fig. EJ). Table also 
includes the values of rj for the 0.33:1 mergers (runs PS- 
PS), which are quite similar to those of the equal mass 
mergers. This shows that the properties of the remnants 
in the projections of the fundamental plane are not very 
sensitive to the mass ratio of the merging galaxies (at 
least for relatively major mergers). 

We can gain an analytical understanding of Fig. Eland 
Table [21 by investigating how a in the R e — M* relation 
varies with the energy and structural parameters rj and / 
in equation 10. Assuming that the stellar mass doubles 
in equal-mass mergers (i.e. M*j = 2M* a very good 
approximation to the simulations) and that the stellar 
effective radius scales as R e oc M" upon merging, equa- 
tion 10 relates the exponent a to rj and /: 

In 2 In 2 ' lj 

We therefore expect a = 1 if a merger is on a parabolic 
orbit with negligible energy transfer (rj = f or b + f t — 0) 
and unchanged structural parameters (// = fi). By con- 
trast, if the merger orbit initially has significant angular 



TABLE 2 

Energetics 



Run" 


r b 


forb 


Tj (= ft + forb) 




fi c 


ff 


PI 








0.225 





367 


0.521 


P2 


2.29 





0.366 





367 


0.493 


P3 


12.5 





0.451 





367 


0.494 


PI 


37.8 





0.531 





367 


0.480 


P5 








0.119 





295 


0.488 


P6 


2.31 





0.317 





295 


0.457 


P7 


12.5 





0.376 





295 


0.463 


P8 


20.5 





0.445 





295 


0.460 


Bl 





0.145 


0.275 





367 


0.513 


B2 


2.29 


0.145 


0.461 





367 


0.497 


B3 


20.8 


0.145 


0.539 





367 


0.488 


B4 


37.8 


0.145 


0.609 





367 


0.456 


B5 


57.1 


0.145 


0.619 





367 


0.433 


B6 


78.4 


0.145 


0.617 





367 


0.423 


a Name of run 



b Pericentric distance of orbit (kpc) 

c For 0.33:1 mergers, runs P5-P8, the value listed is for 
the smaller galaxy (the larger galaxy has the same value 
as in the 1:1 mergers) 



momentum (at distances ^> Ri), the bulges will spiral 
inward and dynamical friction will tend to make the or- 
bits circular at small radii; in this case, rj ss 0.5 by the 
virial theorem and a « 0.58 if ff — fi- In the actual 
simulations, Table [21 shows there is always nonzero en- 
ergy transfer from stars to dark matter (i.e. ft > 0), 
which decreases a for the head-on orbits, whereas the 
merged bulge always has a larger structural parameter 
(i.e. ff/fi > 1), which increases a. Together the vari- 
ations in rj and / explain the behavior and range of a 
(~ 0.7 to 1.3) shown in Fig. 
The fact that ff > fi (by up to ~ 40%) reflects the 
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Fig. 4. — Logarithmic slope a of the stellar mass-size relation 
as a function of the mass ratio of the merging galaxies M1/M2 
calculated from the analytical energy conservation model assum- 
ing ff = Head-on (parabolic) mergers correspond to 77 
while circular orbits have r\ ^ 0.5. The trend of increasing a with 
decreasing orbital angular momentum shown in Fig. 3 for our sim- 
ulations is reproduced in this model, although we do not expect 
the values of a to match exactly due to simplifying assumptions 
made in the model (e.g. ff = fi). The dependence of a on the 
mass ratio is very weak, suggesting that the major merger results of 
Fig. 3 are also applicable to minor mergers. This is also confirmed 
by the similarities between our 1:1 and 0.33:1 simulations. 



varying central dark matter content of the remnants, 
though in general it is non-trivial to quantitatively re- 
late ff/fi to the change in dark matter fraction within 
R e because equation (jHJ) involves a radial integral over 
the stellar density profile and dynamical mass. However, 
from equation |(SJ), we see that for a given distribution of 
stars (fixed R e and p*), increasing the amount of dark 
matter inside the bulge (i.e. increasing Md yn ) increases 
which is why ff > fi in the simulations. 
The analytical energy conservation model given in 
equations J7J and 10 can easily be extended to unequal 
mass mergers. Fig. 0] shows the general dependence of 
a on the mass ratio of the merging bulges for different 
values of 77, assuming for simplicity parabolic orbits and 
// = fi. As argued above, r\ is expected to decrease 
from w 0.5 for orbits with substantial angular momen- 
tum to ~ for nearly parabolic head-on orbits. Corre- 
spondingly, a increases from w 0.58 to 1. This trend of 
increasing a with decreasing orbital circularity (or r pcr i) 
is consistent with that shown in Fig. 3 for the simula- 
tions, although we do not expect the values of a in the 
model and simulations to match exactly due to simplify- 
ing assumptions made in the model (e.g. // = fi). Fig.0] 
also shows that the dependence of a on the mass ratio 
Mi/M 2 is quite weak for 0.01 < M1/M2 < 1.0, indicat- 
ing that the results in Fig. 3 can be extended to very 
minor mergers. This analytic result is consistent with 
the fact that our 0.33:1 mergers and equal mass mergers 
show a similar dependence of a and (3 on merger orbit. 

4. DISCUSSION AND IMPLICATIONS 
4.1. The Fundamental Plane and its Projections 



The results of Section 13.11 show that the fundamental 
plane is preserved during dissipationless mergers for a 
wide range of orbits, a result of virial equilibrium with 
a slight tilt due to an increasing central dark matter 
fraction. On the other hand, the results of Section l3~2l 
show that the projections of the fundamental plane are 
quite sensitive to the amount of energy transferred from 
the bulges to the dark matter haloes during the merger, 
which is determined in large part by the merger orbit. If 
gas-poor mergers help build up the massive end of the 
red sequence, as several current observations and models 
suggest, then our results indicate that the projections of 
the fundamental plane - but not necessarily the plane 
itself - provide a powerful way of investigating the as- 
sembly history of massive elliptical galaxies. 

While a sample of 43 BCGs has been reported to lie 
on the fundamental plane ijOeeerle fc Hoesseil991|) . they 
appear to have different fundamental plane projections 
from those of lower-mass elliptical galaxies. The correla- 
tion between effective radius and luminosity is steeper for 
these BCGs (R e oc L 1 - 25 ) than for all early-type galax- 
ies (R e ex L 6 ) and the BCG velocity dispersions lie 
substantially lower than the extrapol ation of the Faber- 
Jacks on relation to higher luminosity ijOeeerle fc Hoessel 
Il991|) . Other properties of BCGs are notably different 
from those of normal ellipticals a s well, such as their hig h 
incidence of secondary n uclei l| Schneide r et alJ fl983). 
their more prolate shapes l)Rvden et alJll993l). the small 
dispe rsion in their metric luminosities l|Postman fc Laueil 
1995), and the large and relatively constant (out to 
z ~ 0.4) luminosity g ap between BCGs a nd second- 
ranked cluster galaxies ijLoh k, Strau ss 2006). 

The results presented in this paper suggest that these 
differences are all signatures of how massive elliptical 
galaxies were assembled. In particular, dissipationless 
merging of elliptical galaxies provides a natural mecha- 
nism for steepening both the Faber-Jackson and R — L 
relations and preserving the overall fundamental plane, 
provided that merger orbits become preferentially more 
radial for the most massive galaxies (Fig. Ej. 

There is, in fact, support for this possibility in the 
literature. Observationally, it is well-established that 
BCGs are aligned with both their host c lusters and the 
surrounding large scale structure (e.g. iBinggeli Il98"2l 
iFuller et all 119991 iWest & Blakesledl2000|) . Numerical 
simulations sh ow that anisotropic merging is natur al in 
CDM models l)Dubinskil Il99cl IZentner et all I2005D and 
preferentially radial merging along filaments can produce 
many of the observed propert i es of clusters an d BCGs 
iTWesi] IT9941 IWest et a,U IV995t IDubinskil IT99rl . These 
lines of reasoning favour radial mergers during BCG for- 
mation and point to a strong connection between the 
characteristics of the mergers that form BCGs and the 
origin of observed BCG properties (and the differences 
between BCGs and normal ellipticals). In addition, the 
most luminous elli pticals in the Virgo clus ter all seem to 
lie on a filament ( {West &; Blakesleeli2000i). which likely 
means that anisotropic merging is also important for lu- 
minous ellipticals that are not BCGs. 

According to our calculations, projections of the fun- 
damental plane provide independent information about 
the assembly history of massive ellipticals. Analysis of 
a large spectroscopic sample of massive galaxies would 
provide significantly stronger constraints on changes in 
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the fundamental plane projections at the highest masses; 
we expect that the indications of deviations seen in cur- 
rent samples will become more apparent with additional 
data. 

4.2. Supermassive black holes 

Any departure from the canonical Faber-Jackson re- 
lation in massive galaxies would also have significant 
implications for the demography of supermassive black 
holes. In particular, under the assumption that black 
holes coalesce fo llowing a gas-poor merger (wh ich is not 
guaranteed; see iMerritt fc Milosavlievia 120041 for a re- 
view), the black hole mass should increase roughly in 
proportion to the galaxy mass since the amount of mass- 
energy carried aw ay in gravitational waves is likely to be 
quite small (e.g., iBaker et al.1 12004') . We therefore ex- 
pect that dissipationlcss mergers will maintain a linear 
Mbh — M* relation. Any steepening of the L — a c re- 
lation for luminous ellipticals, however, should then be 
reflected in a corresponding stee pening in the M rh — & 
relation for massive black holes. iWvithel l)2006l) has ar- 
gued that just such a steepening is present in the current 
sample of black hole masses, but the statistics are rather 
poor and a larger sample is needed in order to make 
definitive statements. If correct, black hole masses for 
very lumino us galaxies based on the canonical Mrh — 
ct p relation iGebh ardtgl^ al.1 120001: iFerrarese fc Merrittl 
2000; Trcmainc ct al. 2002) would systematically under- 
predict the true black hole mass (see also Lauer et al. 
2006, in prep.). 

4.3. Further discussion 

Several issues can be investigated in more detail to 
test our suggestion that preferentially radial merging is 
responsible for many of the morphological and kinemat- 
ical properties of BCGs. O ne promising possib ility is to 
study the shapes of BCGs. iPorter et alJ (|190 If) find that 
the outer regions of BCGs are generally prolate, wit h the 
inner regions being rounder. iRvden et alJ l) 1993ft find 
that the distribution of BCG ellipticities has a similar 
mean value to that of all ellipticals but with a smaller 
dispersion and a tendency to be slightly more prolate. 
For comparison, we find that the remnants of 1:1 merg- 
ers on low angular momentum orbits are prolate with a 
tendency toward triaxiality: generally a:b:c ~ 1.8:1:1 - 
1.7:1.1:1. By contrast, the runs with significant angular 
momentum have a:b:c ~ 1.7:1.3:1. However, the main 
difference between our equal mass and 0.33:1 mergers is 
that the 0.33:1 remnants are notably rounder than thei r 
1:1 counterparts, an effect also seen bv lVillumsenl Jl982). 
Unequal mass major mergers thus lead to the same over- 
all trends we see in the fundamental plane and its pro- 
jections while producing more spherical remnants. Using 
shapes to quantitatively test the predictions of our mod- 
els will thus require detailed information about the orbits 
and mass ratios of the mergers that formed BCGs. 

The above considerations highlight the need to move 



beyond non-cosmological major merger simulations of 
spherical and isotropic galaxies that have been the fo- 
cus of this paper as well as many earlier papers. In- 
deed, the relaxation and virialization of a cluster gen- 
erally in volves several m assive galaxies with a range of 
masses (.Dubinski 1998). The hierarchical build-up of 
clusters a nd BCGs also likely depends on the clu ster 
mass (e.g. IBrough et allEM lLoh~fc Straiisdl200"3) . In 
addition, the shapes of cluster-mass dark matter haloes 
in dark matter-only cosmological simulations are typi- 
cally quite prolate at formation but then become rounder 
(but still triaxial) a s time goes on l)Hopkins et al]l2005l: 
lAllgood et al.ll2005Tl d ue to a transition fro m filamentary 
to spherical accretion l)Allgood et all2005|) . Clearly, high 
resolution two-component (stars and dark matter) sim- 
ulations in a cosmological context with full formation 
histories are necessary for making detailed predictions 
about the properties of massive ellipticals. We are cur- 
rently working on such simulations. Although we expect 
that our basic conclusions about the fundamental plane 
and its projections will remain unchanged, repeated ma- 
jor and minor merging on realistic cosmological orbits 
will significantly influence the shapes and kinematics of 
the remnant galaxies. 

As one final implication of the non-spherical nature 
of our merger remnants, we note that the measured ve- 
locity dispersions can change by more than 15% with 
viewing angle for a given remnant. For the most mas- 
sive galaxies, this could be significant: a galaxy that has 
cr e w 350 km/s from most viewing angles (typical for a 
BCG) could have er e > 400 km/s in certain rare config- 
urations (we find that the observed velocity dispersion 
is as 15% higher than the mode » 3% of the time for 
the remnants of relatively radial mergers). These ob- 
jects, which are prolate merger remnants viewed down 
the major axis, also have smaller effective radii than most 
viewing angles of rem nants with the same stellar mass. 
iBernardi et al.l l)2005|) have recently reported the detec- 
tion of a population of luminous galaxies with these prop- 
erties. A simple test of our explanation for these galax- 
ies is to look at their projected shapes: our calculations 
predict these galaxies should appear statistically rounder 
than other galaxies of comparable luminosity. This view- 
ing angle effect would also be a source of scatter in the 
Mbh — cr e and Faber-Jackson relations. 
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